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Interacting Chaplygin gas 
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We investigate a kind of interacting Chaplygin gas model in which the Chaplygin gas plays the 
role of dark energy and interacts with cold dark matter particles. We find that there exists a stable 
scaling solution at late times with the Universe evolving into a phase of steady state. Furthermore, 
the effective equation of state of Chaplygin gas may cross the so-called phantom divide w — —1. 
The above results are derived from continuity equations, which means that they are independent of 
any gravity theories. Assuming standard general relativity and a spatially flat FRW metric, we also 
find the deceleration parameter is well consistent with current observations. 
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I. INTRODUCTION 



The existence of dark energy is one of the most signif- 
icant cosmological discoveries over the last century 1]. 
However, although fundamental for our understanding of 
the Universe, its nature remains a completely open ques- 
tion nowadays. Various models of dark energy have been 
proposed, such as a small positive cosmological constant, 
quintessence, k-essence, phantom, holographic dark en- 
ergy, etc., see for recent reviews with fairly complete 
lists of references of different dark energy models. 

Recently the so-called Chaplygin gas, also dubbed 
quartessence, was suggested as a candidate of a unified 
model of dark energy and dark matter 0. The Chaply- 
gin gas is characterized by an exotic equation of state 
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(1) 



where A is a positive constant. The above equation of 
state leads to a density evolution in the form 



Pch 



B 



(2) 



where B is an integration constant. The attractive fea- 
ture of the model is that it naturally unifies both dark 
energy and dark matter. The reason is that, from |(5J, 
the Chaplygin gas behaves as dust-like matter at early 
stage and as a cosmological constant at later stage. Some 
possible motivations for this model from the field theory 
points of view are investigated in PS. The Chaply gm gas 
emerges as an effective fluid associated with c?-branes_ 
and can also be obtained from the Born-Infeld action 
Recently, the original Chaplygin gasmodel was gener- 
alized, with possible observational constraints on these 
generalized models presented in Ref. For example 

in the generalized Chaplygin gas (GCG) approach 
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the equation of state to describe the background fluid is 
generalized to 

A 
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(3) 



and the corresponding evolution of the scale factor is 
given by 



Pch = 



A- 



B 



7,3(1-1-0) 



(4) 



From the above equations, it is clear that when a = \ the 
GCG model recovers the original Chaplygin gas model. 
This approach has been thoroughly investigated for its 
impact on the 0th order cosmology, i.e., the cosmic 
expansion history (quantified by the Hubble parame- 
ter -ff [z]) and corresponding spacetime-geometric observ- 
ables. An interesting range of models was found to be 
consistent with SN la data ■g^, CMB peak locations [3 
and dimensionless coordinate distances to type la super- 
novae ill. There seems to be, however, a flaw in unified 
dark matter (UDM) models that manifests itself only on 
small (Galactic) scales and that has not been revealed 
by the studies involving only background tests. In Ref. 
[Ill, it is found that GCG model produces oscillations 
or exponential blowup of the matter power spectrum in- 
consistent with observations. In fact, from this analysis, 
99.999 % of previously allowed parameter of GCG model 
has been excluded (see, however, 12j). 

Hence we may turn to a model with Chaplygin gas and 
dark matter. Although non-minimal coupling between 
the dark energy and ordinary matter fluids is strongl y re - 
stricted by the experimental tests in the solar system [l^ , 
due to the unknown nature of the dark matter as part of 
the background, it is possible to have non-gravitational 
interactions between the dark energy and the dark mat- 
ter components, without conflict with the experimental 
data. In this paper we investigate some physical prop- 
erties of an interacting Chaplygin gas model. Here, by 
considering an interaction term between the Chaplygin 
gas fluid and dark matter particles similar to those stud- 
ied in the context of quintessence scenarios (15| . we in- 
vestigate dynamical aspects of this interacting Chaplygin 
gas model. 
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Following the more accurate data a more dramatic re- 
sult appears: the recent analysis of the type la super- 
novae data indicate that the time varying dark energy 
gives a better fit than a cosmological constant, and in 
particular, the equation of state parameter w (defined 
as the ratio of pressure to energy density) crosses —1 at 
some low redshift region from above to below 16], where 
w = — 1 is the equation of state for the cosmological con- 
stant. It deserves to note that there are other indepen- 
dent fittings imply the probability that current w < —1 
except for supernovae data The dark energy with 

w < — 1 is called phantom dark energy > for which all 
energy conditions are violated. To obtain w < —1, scalar 
field with a negative kinetic term, may be a simplest 
realization 19]. However, the equation of state of phan- 
tom scalar field is always less than —I and can not cross 
— I. Also it has been shown that the equation of state 
cannot cross —1 in the k-essence model of dark energy 
under some reasonable assumptions 20]. Some dark en- 
ergy models which contain a negative-kinetic scalar field 
and a normal scalar field have been considered in |2ll ] ; in 
these models crossing the border w = —1 can be realized. 
Some different suggestions on this crossing behavior are 
presented in p^ . 

It has been pointed out that Chaplygin gas model can 
be described by a quintessence filed with well-connected 
potential 0. So a model with mutually independent 
Chaplygin gas and dark matter is essentially a special 
quintessence model. The Chaplygin gas , here as dark en- 
ergy, can not cross the phantom divide like quintessence. 
In this paper we shall see that an interaction term can 
realize this crossing naturally. At the same time we ob- 
tain a scaling solution: It may also shed light on the 
coincidence problem. Another interesting result is that 
the scaling solution inevitably leads to the steady state 
Universe 23] , which had been suggested many years ago 
but soon surpassed by expanding Universe, as the final 
state of our model. 

We present our model in details in the next section 
and some observational predictions of this scenario and a 
comparison with recent observational data are also briefly 
discussed. Our conclusions and discussions appear in the 
last section. 



II. THE MODEL 

We consider the original Chaplygin gas, whose pressure 
and energy density satisfy the relation, pch = —A/pch- 
By assuming the cosmological principle the continuity 
equations are written as 



Pch + SHjchp 



'ch 



-r. 



and 



Pdm + SHjdmPdn 



(5) 



(6) 



where the subscript dm denotes dark matter, H is the 
Hubble parameter and 7 is defined as 



7 



1 



P 
P 



1 



■ w, 



(7) 



in which w is the parameter of the state of equation, 
and = 1 throughout the evolution of the Universe, 
whereas "fch is a variable. 

r is the interaction term between Chaplygin gas and 
dark matter. Since there does not exist any microphys- 
ical hint on the possible nature of a coupling between 
dark matter and Chaplygin gas (as dark energy), the in- 
teraction terms between dark energy and dark matter are 
rather arbitrary in literatures |24| . Here we try to present 
a possible origin from fundamental field theory for T (see 
|25l ] for a thermodynamic discussion on F). 

Whereas we are still lacking a complete formulation of 
unified theory of all interactions (including gravity, elec- 
troweak and strong) , there at present is at least one very 
hopeful candidate, string/M theory. Although the recent 
developments in string theory, assisted by the discovery 
of the power of duality, have greatly improved our un- 
derstanding of it, the theory is still not known in a way 
that would enable us to ask the questions about space- 
time in a general manner, say nothing of the properties 
of realistic particles. Instead, we have to either resort 
to the effective action approach which takes into account 
stringy phenomena in perturbation theory, or we could 
study some special classes of string solutions which can be 
formulated in the non-perturbative regime. But the lat- 
ter approach is available only for some special solutions, 
most notably the BPS states or nearly BPS states in the 
string spectrum: They seems to have no relation to our 
realistic Universe. Especially, there still does not exist 
a non-perturbative formulation of generic cosmological 
solutions in string theory. Hence nearly all the investi- 
gations of realistic string cosmologies have been carried 
out essentially in the effective action range [2^. Note 
that the departure of string-theoretic solutions away from 
general relativity is induced by the presence of additional 
degrees of freedom which emerge in the massless string 
spectrum. These fields, including the scalar dilaton field, 
the torsion tensor field, and others, couple to each other 
and to gravity non-minimally, and can influence the dy- 
namics significantly. Thus such an effective low energy 
string theory deserve research to solve the dark energy 
problem. There a special class of scalar-tensor theories 
of gravity is considered to avoid singularities in cosmolo- 
gies in ^27]. The action is written below. 
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where G is the Newton gravitational constant, is a 
scalar field, L,„ denotes Lagrangian of matter , ^ repre- 
sents different matter degrees of matter fields, q guaran- 
tees the coupling strength between the matter fields and 
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the dilaton. With action (jSJ , the interaction term can be 
written as follow l27ll. 



Hpr 



d\nq' 
d In a 



(9) 



Here we introduce new variable q{a)' = (/(a)^'^'""^^-'/^, 
where a is the scale factor in standard FRW metric. By 
assuming 



(10) 



where pm and are the densities of matter and the scalar 
field respectively, one arrive at the interaction term, 



r = 3Hc{p,n + P?), 



(11) 



which is just the coupling form studied in contexts 
of quintessence and phantom dark energy models |l5j . 
Moreover the Chaplygin gas can be view as a scalar 
field with proper potential in cosmological models 
So it may be reasonable to phenomenologically introduce 
such an interaction term between Chaplygin gas and dark 
matter. 

Now return to the equation set ^ and 0. Set 
X ^ - ln(l+z), r = 3Hc{pch+Pdm), u = {iHl)-^ pdm, 
V = {-iHiy^i^pch, A' = A{ml)-'^K^, where Hq takes 
the value of present Hubble parameter, k, is the New- 
ton gravitational constant and c is a constant without 
dimension. Eqs. (jSJ and 10 reduce to 



du 
dx 



-iu + 3c(7i + v), 



(12) 



dv 

= A' /v)-ic{u + v). (13) 

dx 

We note that the variable time does not appear in the 
dynamical system H12|) and H13() because time has been 
completely replaced by redshift x — —ln{l + z). The crit- 
ical points of dynamical system (|12|) and (|13|l are given 

by 



du dv ^ 
dx dx 

The solution of the above equation is 



Ur = 



-Vr 



1 - C 

,l = {l-c)A'. 



(14) 

(15) 
(16) 



We see the final state of the model contains both Chap- 
lygin gas and dark matter of constant densities if the 
singularity is stationary. The final state contents perfect 
cosmological principle: the Universe is homogeneous and 
isotropic in space, as well as constant in time. Physically 
r in (0) plays the role of matter creation term C in the 
theory of steady state universe at the future time-like in- 
finity. Recall that c is the coupling constant, may be pos- 
itive or negative, corresponds the energy to transfer from 



Chaplygin gas to dark matter or reversely. A' must be a 
positive constant, which denotes the final energy density 
if c is fixed. Also we can derive an interesting and simple 
relation between the static energy density ratio 



l+Ts 



where 



1 . Pdm 

Z — Pch 



(17) 



(18) 



To investigate the properties of the dynamical system in 
the neighbourhood of the singularities, impose a pertur- 
bation to the critical points. 



d{5u) 
dx 



-iSu + iciSu + Sv), (19) 



'^!^^-^5v + ^6v)-ic{5u + 5v). (20) 
dx 

The eigen equation of the above linear dynamical system 
(Su, Sv) reads 

(A/3)2 + (2+-^)A/3 + 2-2c2 = 0, (21) 
1 — c 



whose discriminant is 



A = [(1 - c)* + (3/2 - c)2]/(l - cf > 0. 



(22) 



Therefore both of the two roots of eigen equation 1)21(1 
are real, consequently centre and focus singularities can 
not appear. Furthermore only G (0,cx3), such that 
c G (0, 1), makes physical sense. Under this condition it 
is easy to show that both the two roots of (|21|) are neg- 
ative. Hence the two singularities are stationary. How- 
ever it is only the property of the linearized system H12() 
and ((T^ , or the property of orbits of the neighbourhoods 
of the singularities, while global Poincare-Hopf theorem 
requires that the total index of the singularities equals 
the Euler number of the phase space for the non-linear 
system (|19|l and (|20ll . So there exists other singularity 
except for the two nodes. In fact it is a non-stationary 
saddle point at u = 0, v = with index —1. This sin- 
gularity has been omitted in solving equations (|12|) and 
(|13|l . The total index of the three singularities is 1, which 
equals the Euler number of the phase space of this plane 
dynamical system. Hence there is no other singularities 
in this system. From these discussions we conclude that 
the global outline of the orbits of this non-linear dynam- 
ical system (|12|l and (|13|) is similar to the electric flux- 
lines of two negative point charges. Here we plot figures 
1 and 2 to show the properties of evolution of the Uni- 
verse controlled by the dynamical system H12I) and H13|l . 
As an example we set c — 0.2, A' = 0.9 in all the figures 
except figure 0] 

Further, to compare with observation data we need the 
explicit forms of u{x) and v(x), especially v{x), since we 
have set 7^™ = 1 but jch is not a constant. We need 



FIG. 1: The plane v versus u. (a) We consider the evolution of the universe from redshift z = — 1. The initial condition 
is taken as u = 0, v — 400; u = 50, v = 350; u — 100, v — 300; u = 120, v = 280 on the four orbits, from the left to 
the right, respectively. It is clear that there is a stationary node, which attracts most orbits in the first quadrant. At the 
same time the orbits around the neibourhood of the singularity is not shown clearly, (b) Orbit distributions around the node 

Uc = VcC/(l - c), Vc = V'(l - C)A'. 



the properties of jch in our model, which is contained in 
v{x), to compare with observations. EUminate u{x) by 
using ifT^ and ((T^ we derive 



1 d-'^'" r-, /, A/ , 9\ , ^ dv 

+ + + )/^]:r 

6c ax'^ ax 



3(1 



dv 
dx 



3(i; - A'/v) 



/(3c)l=0, (23) 



which has no analytic solution. We show some numerical 
solutions in figure |31 We find that for proper region of 
parameter spaces, the effective equation of state of Chap- 
lygin gas crosses the phantom divide successfully. We 
have no analytical result on this crossing phenomena yet 
in the present stage. 

Up to now all of our results do not depend on Einstein 
field equation. They only depend on the most sound prin- 
ciple in physics, that is, the continuity principle, or the 
energy conservation law. Different gravity theories, such 
as standard general relativity, brane-induced gravity, l/R 
gravity, and Lovelock gravity or cosmology of mod- 
ified Friedmann equation such as Cardassian cosmology 
|29j . correspond to different constraints imposed on our 
previous discussions. Our improvements show how far we 
can reach without information of dynamical evolution of 
the Universe. 

The most significant parameters from the viewpoint of 
observations is the deceleration parameter g, which car- 
ries the total effects of cosmic fluids. From now on we 
introduce the Friedmann equation of the standard gen- 
eral relativity. As a simple case we study the evolution 
of g in a spatially flat frame. So q reads 



n-2 



1 fu + v-^A'/v'^ 

2 V u + v 



(24) 



ter V should satisfy 



u{0)+v[id) = 1. 



(25) 



And then Friedmann equation ensures the spatial flatness 
in the whole history of the Universe. Before analyzing the 
evolution of q with redshift, we first study its asymptotic 
behaviors. When z oo^ q must go to 1/2 because both 
Chaplygin gas and dark matter behave like dust , while 
when z — !■ — 1 g is determined by 



lim q 



iA'/v^ 



(26) 



One can finds the parameters c = 0.2, A' — 0.9 are dif- 
ficult to content the previous constraint Friedmann con- 
straint (|25|l . Here we carefully choose a new set of pa- 
rameter which satisfies Friedmann constraint H25|l . say, 
A' = 0.4, c = 0.06. Therefore we obtain 



lim q 

2—^ — 1 



-1.95, 



(27) 



and density of Chaplygin gas u and density of dark mat- 



by using l(TK|l and (fT?)|l . Then we plot figure 0] to clearly 
display the evolution of q. One can check u{fi) = 
0.25, w(0) = 0.75; u(0) = 0.28, t;(0) = 0.72; u(0) = 
0.3, w(0) = 0.7, respectively on the curves v{—2) = 
273; w(-2) = 250; t;(-2) = 233. One may find an inter- 
esting property of the deceleration parameter displayed 
in figure ^ The bigger the proportion of the dark en- 
ergy, the smaller the absolute value of the deceleration 
parameter. The reason roots in the extraordinary state 
of Chaplygin gas in which the pressure p^h is in- 
versely proportional to the energy density pch- 

Also we note that maybe an FRW Universe with non- 
zero spatial curvature fits deceleration parameter better 
than spatially flat FRW Universe. This point deservers 
to research further. 
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FIG. 2: The plane v versus u. (a) To show the global prop- 
erties of dynamical system ^Vl\ and we have to include 
some "unphysical " initial conditions, such as u = —100, v — 
—300, except for physical initial conditions which have been 
shown in figure^ (b) Orbits distributions around the nodes. 
The two nodes Uc = Vcc/{1 — c). 



= ^(1 - c)A' and 



Uc = Vcc/{1 — c), Vc — ~ c)A' keep reflection sym- 

metry about the original point. Just as we have analyzed, we 
see that the orbits of this dynamical system are similar to the 
electric fluxlines of two negative point charges. 
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FIG. 3: V versus x. The initial condition of the three curves is 
u{-2) = 0, v{-2) = 400; u{-2) = 50, v{-2) = 350; u{-2) = 
100, v{ — 2) — 300, respectively. Obviously the energy density 
of Chaplygin gas rolls down and then climbs up in some low 
redshift region. So the Chaplygin gas dark energy can cross 
the phantom divide ui = — 1 in a fitting where the dark energy 
is treated as an independent component to dark matter. 




FIG. 4: q versus x. The initial condition of the three curves 
is -^(-2) = 0, v{-2) = 273; m(-2) = 15, v{-2) = 250; 
u{ — 2) = 25, v{—2) = 233, respectively. Evidently the decel- 
eration parameter q of Chaplygin gas rolls down and crosses 
q = in some low redshift region. The transition from decel- 
eration phase to acceleration phase occurs ai z — 0.18; z = 
0.21; z = 0.23 to the curves u{-2) = 0, v{~2) = 273; 
u{-2) = 15, v{-2) = 250; ■u(-2) = 25, v{~2) = 233, re- 
spectively. One finds —q ~ 0.5 ~ 0.6 at z = 0, which is well 
consistent with observations. 



III. CONCLUSION AND DISCUSSION 

We present a phase-space analysis of the evolution for a 
FRW universe driven by an interacting mixture of dark 
matter and Chaplygin gas. In the absence of interac- 
tion, there exist no scahng solutions because the state 
of equation of Chaplygin gas decreases with scale fac- 
tor while state of equation of dark matter keeps a con- 
stant. Hence we study the existence and stability of the 
cosmological scaling solutions with interaction between 
Chaplygin gas and dark matter. For the interaction term 
r = 3cH{pch + Pdm), inspired by low energy effective 
string theory and scalar-tensor gravity theory, we find 
stationary scaling solutions for reasonable initial condi- 
tions. Further more this approach leads to an impressive 
result: The equation of state of Chaplygin gas traverses 
the phantom divide it; = — 1, which is favored by recent 
fittings with current type la supernovae data. At the 
same time, different from phantom model, there is no 
singularity in the future in our model. On the contrary 
the final state of our model is steady state Universe, in 
which Chaplygin gas ensures the continuous production 
of matter (dust) through the interaction term. The an- 
alytical researches on this crossing behavior deserve to 
investigate in the future work. We also calculate the de- 
celeration parameter of this model in frame of spatially 
flat FRW cosmology. The result is consistent with the 
observation data. 

At last, as a phenomenological model, we must con- 
strain it by observational data. We see that this model 
is so simple that it is fully parameterized by two param- 
eters: c determines the final ratio of Chaplygin gas and 
dark matter, and A governs the final total energy den- 
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sity of the Universe. We may consider the observational 
constraints on the parameter space arising from observa- 
tion data from different observations, such as supernovae, 
CMB, X-rays, gravitational lensing effects or combina- 
tion of these data in future work. 
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